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DISCLAIMEB 



Although each program has been tested by its contributor, no 
warranty, express or Implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USEFS Group, in con- 
nection therewith. 



1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 



(fill out in typewriter or pencil, do not use ink) 
Program No. _______ Date 

Program Name: 



1. Does the abstract adequately describe what the program is and what Yes No 

it does? " ~™ 
Comment 

2. Does the program do what the abstract says ? Yes No 

Commen t "™ 

3. Is the Description clear, understandable, and adequate? Yes No 

Commen t • ™ A " 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No 

Commen t " m ' 

Are the Sense Switch options adequately described (if applicable)? Ye s No 

Are the mnemonic labels identified or sufficiently understandable? Yes No 

Commen t ~™ ~ 

5. Does the source program compile satisfactorily (if applicable)? Ye s No 

Commen t 

6. Does the object program run satisfactorily? Yes No 

Commen t " "™ 

7. Number of test cases run , Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No 

Comment 

8. Does the Program Meet the minimal standards of the 1620 Users 

Group? Yes No 

Comment 

9. Were all necessary parts of the program received? Yes No 

Commen t 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 



Mr. Richard L. Pratt' Company 
Data Corporation 

7500 Old Xenia Pike Address 
Dayton, Ohio 45432 



User Group Code 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATION'S PROGRAM 
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PROGRAM No. lit AND 70U| - NEW ENGLAND ELECTRIC SYSTEM 
Ti.'.Ie i Simultaneous Equation Solution and Matrix Inversion with Complex Values 
Written ty s G. S. Haralampu, New England Electric S 70 tat 

Scope s The program has been written for the solution of a maximum of fifteen 

equations with fifteen unknowns for the basic 1620, with the coefficients, 
unknowns, and constants being coaplex numbers. The matrix inversion aethod 
is used for the solution of the sianltaneous equations* (Rotes k 2$ x 2$ 
matrix is allowed on the 7010 computer). 

Machine Configuration s 20K 1620 with a 1622 card reader and punch, or 
80K 7010 computer 

Reference ; "Numerical Matheamtical Analysis" by Scarborough, Fourth Edition, Page 5U0. 

Speed s A It x It matrix requires 20 seconds for solution. A 16 x 16 matrix requires 12 
minutes for solution, on the 1620. The speed on the 7010 is approximately 
12 times faster than the 1620. 

Analysis s Assume the following k x k matrix is to be inverted: 



All the Z's are complex numbers . 



The following aethod was used: 

A. Place the first colwan of a unit matrix, equal to the siie of 
the main matrix, to the right of the main matrix. This la 
shown belowt 

Main Matrix First Column of the Unit Matrix 



hi 


212 


Z 13 


hh 


Z 21 


Z 22 


Z 23 


Z 2h 


Si 


Z 32 


Z 33 


Z * 


Z Ul 


Z U2 


Z U3 





hi 


h.2 


^3 


hk 


1 


hi 


"hi 


z 23 







Z 31 


Z 32 


Z 33 


Z 3U 





Z U1 


Z U2 


Z U3 


Z Uk 






B. By mathematical manipulation, obeying all the rules of matrix 
handling of rows and columns, obtain the first column of the unit 
matrix at the left of the main matrix. This will bes 



First Column of 
the Unit Matrix 



Intermediate Values 
of the Main Matrix 



1 


*12 


*13 


IlU 


*15 





*22 


*23 


x 2 U 


Y 25 





*32 


*33 


*3U 


Y 35 





*U2 


*U3 







Shift the Intermediate values of the matrix to the left by one column 
and place the second column of the unit matrix to the right. 



Main Matrix Second Column of the Unit Matri; . 



Ill 


*12 


113 


111, 





Y21 


T22 


*23 


*21» 


1 


*31 


*32 


*33 







Da 


ru2 


T U3 








Repeat "B" until we haves 



Second Column of Second Intermediate Values 
the Unit Matrix of the Main Matrix 






*12 


*13 


Y lU 


Y l5 


1 


I 22 


*23 


Y 2l» 


Y25 





*32 


*33 


*3U 


Y35 





Tjj 2 


*1*3 




Y U5 



Repeat parts B and C until the matrix manipulation has been done for ti-.r 
four columns of the unit matrix. Then the values left in position 



*12 


*13 


111* 


H5 


T22 


*23 


Y2U 


*25 


*32 


*33 


T3U 


T35 


%2 


*U3 


Yy, 


*U5 



compose the inverted matrix. 



Once the matrix is inverted, the unknowns can be found by matrix multi- 
plication. Assume the given equations to be solved are in the form of 

[~Z~| x [i] = [e] Equation 1 

where each bracketed symbol represents a matrix. Solving for the I's 
we have: 

(V) = [y] x [jj Equation 2 

where the Y matrix is the inverse of the Z matrix. 

If the equations to be solved are as shown below, 

Zll 1 ! ♦ z 12*2 * z 13 x 3 + Z lU T ii " E l 
z 21 x l * Z 22 J 2 * z 23 x 3 * z 2l» I l» " E 2 
z 31 x l * z 32 x 2 * Z 33 J 3 + Z 3U I U " E 3 
^l 1 ! + Z U2 I 2 * Z U3 I 3 * Wli " E U 
then Equation 1 expanded is as follows: 



z ll 




h2 


Z 13 


hk 




h 




VI 


Z 21 




Z 22 


223 


Z 2U 


x 


h 




*| 


Z 31 




Z 32 


Z 33 


Z 3U 




h 




E3 j 


Z U1 




Z U2 


Z U3 


Z UU 




_\ 






i Equation 2 is: 




















*12 


*13 








El 


12 




121 


Y22 


*23 






x 


E 2 


13 




^31 


v 32 


^33 




Y 3U 




E 3 






Y Ul 


\2 


Y U3 




Y W 




_ E U_ 



Upon completing the matrix multiplation, the unknown I's are: 
l x - Y u Ei ♦ Y 12 E 2 ♦ Y13E3 * Yii,E U 

1 2 - Y 21 E 1 ♦ Y 22 E 2 ♦ Y23E3 ♦ Y 2U E U 

1 3 - Y 31 E 1 ♦ Y 32 E 2 ♦ Y 33 E 3 ♦ Y 3U E 1 , 

\ ' Y U1 E 1 + Y l42 E 2 * Y 1*3 E 3 + Vh 



USERS' GUIDE 
INPUT DATA FOR THE 1620 OR THE 7010 

First Card : Title of job with a minus sign in column 72. This will skip to a new 
page on the U07 far each separate solution. 

Second Card : (Right justified in the first five coliams). The number of equations 
to be solved. (Fixed point) 

Subsequent Cards : (Floating point) The real part of the coefficient is entered in the 
first ten columns, and the imaginary part in the next ten columns. 
For example, for a three equation solution, the input ist 



z ll 


z 12 


^3 


El 










Z 21 


Z 22 


Z 23 


E 2 


Z 31 


Z 32 


Z 33 


E 3 



and 12 cards are required. 

Last Card: (Fixed point; right justified in the first five colxmma) 

a) The value 1 if another case is to follow, and a 9999 if this was the 
last case* 

OUT POT 

See examples in the Appendix. 

LIMITATIONS 

1) Largest matrix handled is a 15 x 15 for the 1620, and 25 x 25 for the 7010. 

2) No seroes allowed in the main diagonal* 

3) If a larger than a 15 x 15 on the 1620 (or a 25 x 25 on the 7010) is attempted 

to be read in, the computer will not solve that case but will print message saying 
so, and then proceed to the next case* 



- 5 - 

FOR THE 7010 ONLY 

Place the following card ahead of the input data: 

(Columns 6-10) M0N$| 

(Columns 16-26) EIEQ HARA,MJB 

A Matrix Inversion Subroutine has also bean wHttjm 7mn 
To make use of this subroutine do the following: * 7 ° 10, 

A. In your win program include the statements: 

1. DIMMSION R(25,26), X (25,26) 

2. COMMON R,X,K 

3. Read K or set K to equal size of matrix to be solved. 

1*. CALL MATINV(R,X,K) 

B * ^JfV* °? th * lnpUt ' inclwto the following card before the »M0II$1 
EXEQ " card of your program: 

(Columns 16-26) CALL MATIN V 



- 6 - 

OPERATING INSTRUCTIONS 

Clear Memory 

Load Object program in computer with the input for all the cases to be solved. 

When solution is completed for all the cases, the 1620 will PAUSE, and the 
7010 will EXIT. 
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INDEX TO THE APPENDIX 

Appendix I - Input and output fron the 1620 

Appendix II - Fortran listing for the 1620 

Appendix III- Output fro* the 7010 

Appendix IV - Fortran listing for the 7010 

Appendix V - Fortran listing for the aatrix inreraion subroutine for the 7010 

Appendix VI - Fortran prograa liatiag solring siamltaneoua equations, using the 
■etrix inreraion subroutine. This mm written to test the MATIHV 
subroutine. 



INPUT DATA 



APPENDIX I 



THIS 


IS A 4X4 MATRIX 


4 




0.071 


.363 


-8.023 


-4.031 


e. 148 


4.856 


8.646 


6.323 


1 .000 


0.0 


.042 


.236 


3.217 


7.501 


-3.31 1 


-7.813 


• 665 


2.632 


1 .000 


0.0 


-1 .4S2 


-5,067 


-9.569 


-9.686 


8.050 


4.390 


10.167 


12.109 


1.00 


6.0 


.097 


• 564 



-B.lOl -4,564 
8.105 4.718 



8.715 7.022 
1_.00 0.0 
1 

THIS IS * 3X3 MATRIX 
3 

.0582 .168 



•0 -.00 19* 

,0 -.000465 

,058 " -.00502 

0,0 -,00196 

To55I . 1 68 

0,0 -,00196 

'"-•1102 ,6000687 

0»0 -.000465 



0,0 -,00196 
.0582 .168 



-,1102 .0000687 
9999 



APPENDIX H 



57070 
57070 
67070 
67070 
57070 
57070 
57070 
57092 

57130 
57252 
57384 
57504 
57588 
57666 
5777 1 * 
67896 
57936 
579^8 
57948 
57972 
58040 
08052 
58064 
58256 
58376 
58376 
68376 
68400 
58436 
58448 
58460 
58544 
58664 
58748 
68772 
58852 
58932 
59112 
59160 
59208 
09220 
59436 
59664 
59748 
59868 
59892 
59916 
59984 
T0020 
T0044 
T0052 
T0088 
Toil 2 
T0124 



1000 
1001 
1002 
1003 
1004 
1005 
1006 
1007 
1008 
1009 
55 



103 



30 
50 



48 



31 
34 



35 

38 

37 
39 



PROGRAM NUMBER 14**N0V 1,1964 

SIMULTANEOUS EQUATION SOLUTION USING MATRIX INVERSION METHOD 
COMPLEX NUMBERS ARE USED 

WRITTEN BY G .S .HARALAMPU - NEW ENGLAND ELECTRIC SYSTEM 
SIZE IS LIMITED TO A 15X15 MATRIX 

DIMENSION R(15,16).X(15.16),ER(15),EI(15).TR(15).TI (15) 

FORMAT ( 15) 

F0RMAT(2F10.4,2I5) 

FORMAT (/46H THE INVERTED MATRIX ELEMENTS ARE AS FOLLOWS) 

FORMAT (// 46H THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW ) 

FORMAT (48HMATR IX LARGER THAN 15X1 5. NEXT CASE TO BE SOLVED.) 

FORMAT (30H REAL I MAG ROW COL) 

FORMAT (27H REAL I MAG ELEMENT) 

FORMAT (42H AN ELEMENT ON THE MAIN DIAGONAL IS ZERO) 

FORMAT (49H ) 

FORMAT(36X,35X,1H-) 

READ 1008 

READ IN THE SIZE OF THE MATRIX 

READ 1000.K 

IF(K-16)103,101 ,101 

DO 50 1=1 ,K 

DO 30 J m 1, K 

READ 1001, R(I,J),X(I,J) 

READ 1001, ER(I),EI(I) 

READ IN THE LAST CODE CARD. 1 FOR MORE CASES TO FOLLOW, AND 9999 

FOR THE LAST CASE 

READ 1000.NCOOE 

KX - K + 1 

DO 33 I - 1, K 

DO 48 L-1.K 

R(L, KX) - 0.0 

X(L, KX) - 0.0 

R(I,KX)«1.0 

IY - 1 

IF(R(I,1))34,31.34 

IF(X(I,1))34,100,34 

T1-R(I,1)*R(U1)+X(I,1)*X(I ,1) 

T1R-R0 J) 

T1X-X(I,1) 

DO 35 J-1.KX 

WR«(R(I,J)*T1R+X(I ,J)*T1X)/T1 
WU(X( I ,J)*T1R-R( I ,J)*T1X)/T1 
R(l , J) - WR 
X( I , J) - Wl 

IX -0 

IY - 2 

IF(I-K)37.38,37 

MX- 1-1 

MY-1 

GO TO 39 
MY m I + 1 
MX - K 

DO 40 L - MY, MX 



IX 



IX +1 



To 160 

T0208 
T0256 
T0268 

T0484 
T0688 
T0844 
T1072 
T1096 
Til 64 
T1244 
T1256 
T1268 
T1304 
T1448 
T1664 
T1700 
T1712 
T1760 
T1808 
T1820 
T2084 
T2336 
T2420 
T2576 
T2588 
T2600 
T2612 
T2624 
T2636 
T2648 
T2900 
T2912 
T2924 
T2936 
T3068 
13 136 
T3148 
T3160 
T3172 
13180 

13192 
3204 
3216 
T3288 
13360 
T3384 
T3392 

T9979 
T9949 
T9909 
T9869 



40 



53 



49 



201 
202 
102 

100 

101 



2001 
2000 



JO 

T1 - R(L, 1) 
T2 = X(L. 1) 
DO 40 J«1,KX 

WR - R(l , J)*T1 - X( I , J)*T2 
Wl - R(l, J)*T2 + X(l, J)*T1 
R(L, J) = R(L, J) - WR 
X(L,J)-X(L,J)-WI 
IY m 3 

IF( 1-1)45,44,45 

IF((K-1)-IX)38,44,38 

DO 46 L - 1 , K 

DO 46 J - 1, K 

NU - J + 1 

R(L,J)-R(L,NU) 

X(L, J) - X(L, NU) 

CONTINUE 

DO 53 1-1, K 

TR(I )-0.0 

TKD-O.O 

DO 53 J»1,K 

T1-R(I,J)*ER(J)-X(I ,J)*EI(J) 

T2-X(I,J)*ER(J)+R(I,J)*EI(J) 

TR(I)-TR(I)+T1 

TI(I)-TI(I)+T2 

PUNCH 1009 

PUNCH 1008 

PUNCH 1002 

PUNCH 1005 

00 49 l»1 ,K 

00 49 J - 1, K 

PUNCH 1001, R(I,J),X(I,J),I,J 
PUNCH 1003- 
PUNCH 1006 
DO 201 1-1 ,K 

PUNCH 1001, TR(1 ),TI(I),I 

IF(NC0DE-1)102,55,102 

PAUSE 

PAUSE 

PRINT 1007 
GO TO 202 
PRINT 1004 
00 2000 1-1, K 
00 2001 J-1.K 
READ- 1001 ,W,Z 
READ 100t,W,Z 
READ 1 000, N CODE 
GO TO 202 
END 



COS 
ATANF 
LOGF 
SIGNF 



19939 EXP 
19899 SQRT 
T9859 ABS 



T9999 SIN 

59969 COSF 
9929 EXPF 
9889 SQRTF 
9849 ABSF 



[9989 SINF 
19959 ATAN 
[9919 LOG 
[9879 SIGN 
F9839 R T7449 



© 



a 



T*739 TR 
T**29 T001 
T*389 T005 
T009 
TU309 OOO 
1*269 I 
T*229 KX 

t*1*9 503* 
TU109 001 

Itofi Sod* 

T3989 «X 
T39*9 J2 
T3909 50*6 
T3869 6201 
T3829 2*001 



17*39 X 
T*599 



T50*9 



LOAD SUBROUTINES 



T**19 T002 
1*379 1006 
T*339 5055 
1*299 5l03 
1*259 §010 
lklX9 BOOT 
1*179 6000000055 
f*U9 003 1 
TM>99 TIR 
1*059 §02 v 
t*0t9 OOOl 
13979 MV . 
13939 0005 
T3899 «U 
T3859 5202 
T3819 W 



T5039 
T*S89 
T**09 
1*369 
T*329 
J*289 
1*2*9 
T*209 
t*U9 
t*l29 
1*089 
f*0*9 
1*009 
13949 
13929 
13889 
138*9 
13809 



ER T*899 
Tl JkkkS 

T003 
I007 

K 

BlOl 

J 

5033 ^ 

TOOOOOOOSl 

0100 

XIX 
Wl 

9037 
8039 



S3! 

S102 



T*889 
T**39 
T*399 
1*359 
1*319 
1*279 
1*239 
1*199 



El T*7*9 

Tooo 

Too* 

T008 

6otS 

6oso 

NCOOE 
50*8 




A n PENra:x in 



OUTPUT 



THIS A 4X4 MATRIX 



THE INVERT 
REAL 

-.1862 

-.0085 

-.01*3 
.2154 
.2298 
.029? 

-.0062 

-.1673 
.2177 

-.0080 

-.0056 

-.1386 
.0327 
•0197 
.0054 
.0175 



ED MATRIX ELEMEN 



I MAG 
1.0951 
.0151 
.1282 
-1.2256 
-1.2040 
-.0935 
.0331 
1.1311 
-1.1513 
-.0049 
.0339 
1.0672 
-.0725 
-.0779 
-.0046 
.0319 



ROW 
1 
1 
1 
1 
2 
2 
2 
2 
3 
3 
3 
3 
4 
4 
4 
4 



CO 



ARE AS FOLLOWS 



THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW 
REAL I MAG ELEMENT 



-.0117 


.0129 


I 


.0858 


-.1339 


2 


.0684 


-.0351 


3 


•0794 


-.1231 


4 



A^rrax in 

OUTPUT 



THIS IS A 3X3 MATRIX 

THE INVERTED MATRIX ELEMENTS ARE AS FOLLOWS 



REAL 


I MAC 


ROW 


COL 


1.6*17 


-5.3149 


1 


1 


.0385 


-.0488 


1 


2 


.0096 


-.0119 


1 


3 


.0385 


-.0488 


2 


1 


1.8*22 


-5.3153 


2 


2 


.0385 


-.0488 


2 


3 


.0096 


-.0119 


3 


1 


.0385 


-.0488 


3 


2 


1.8417 


-5.3149 


3 


3 



THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW 
REAL I MAG ELEMENT 

.0748 -.3108 1 

-.2049 .5882 2 

-.2063 .5904 3 



C PROGRAM NUMBER 14«*MAY 1.1964 

C SIMULTANEOUS EQUATION SOLUTION USING MATRIX INVERSION METHOD 

C COMPLEX NUMBERS ARE USEO 

C WRITTEN BY G. S.HARALAMPU - NEW ENGLAND ELECTRIC SYSTEM 

C SIZE IS LIMITED TO A 25X25 MATRIX 

DIMENSION R < 25 • 26 ) » X ( 25 * 26 ) • ER < 25 ) • E I < 25 ) » TR < 25 ) ♦ T I ( 25 ) 

1000 FORMAT ( 15) 

1001 FORMAT (2FI 0.4. 215 ) 

1002 FORMAT (/46H THE INVERTED MATRIX ELEMENTS ARE AS FOLLOWS) 

1003 FORMAT < ///46H THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW ) 

1004 FORMAT (67H THIS MATRIX IS GREATER THAN 25X25. THE NEXT CASE WILL 
1 BE SOLVED. ) 

1005 FORMAT OOH REAL I MAG ROW COL) 

1006 FORMAT (27H REAL I MAG ELEMENT) 

1007 FORMAT (42H AN ELEMENT ON THE MAIN DIAGONAL IS ZERO ) 

1008 FORMAT (IX. 52H > 
3000 FORMAT < 1H 1 ) 

55 READ (1.1 008 ) 
C READ IN THE SIZE OF THE MATRIX 

READ (1.1 000 ) K 
IF(K.GT.2S)GO TO 101 
DO 50 I«1«K 
DO 30 J » 1. K 
30 READ(l.IOOl) R(I.J).X(I»J) 
50 RE AO (1.1001) ER(I).EKI) 
C REAO IN THE LAST CODE CARD. 1 FOR MORE CASES TO FOLLOW. AND 9999 

C FOR THE LAST CASE 

READ (1.1 000 ) NCODE 
KX » K + 1 
DO 33 I « I» K 
DO 48 L*l «K 
R(L. KX) » 0.0 
48 X(L. KX) « 0.0 
R( I «KX>*1 .0 
IV » 1 

IF(R(I.l ).NE.0.0)GO TO 34 
IF(X( 1.1 ).EQ.O.0)GO TO 100 

34 Tl«R(I.l)«R(t.l>*X(I.l )*X(I.l ) 
T1R>R(I. 1 ) 

T1X«X<I.1 ) 
DO 39 J*l .KX 

WR*(R< I . J)«T1R+X( I . J)»T1X)/T1 
WI««X( I . J)«T1R-R( I«J)«T1X)/T1 
R ( I • J > » WR 

35 X(l, J) • «! 
IX «0 

IY ■ 2 

IF(I.NE.K)GO TO 37 

38 MX« I - 1 
MY«1 

GO TO 39 
37 MY ■ I ♦ 1 
MX - K 

39 DO 40 L * MY. MX 
IX - IX +1 

Tl ■ R(L. 1 > 
T2 ■ X(L. 1 ) 



o 



DO 40 JM.KX 

WR » R<I. J)*TI - X<I. J)*T2 
MI s R(I. J)*T2 + XU. J)*T1 
R (L • J) * R(Li J) - WR 
40 X(L. JI«X(L. J)-WI 
IY * 3 

IF< I.EQ.l )GO TO 44 

IFMK-1 ).NE» IX )60 TO 38 
44 DO 46 L • 1 » K 

DO 46 J > 1 < K 

NU « J ♦ 1 

R(L*J)>R<L.NU) 
46 X(L. J> - X(L, NU) 
33 CONTINUE 

DO 93 I-ItK 

TR( I >«0.0 

TI ( I )"0.0 

DO 53 J-l.K 

T 1 "R ( I »J)«ER< J)-XU «J)»EI (J) 

T2»X« I.J) *ER < J ) +R ( I . J)*EI (J) 

TR< I )«TR< I )+Tl 
S3 TI < I »"TI ( I »+T2 

WRITEO.3000) 

WRITEO.1008) 

WRITE<3.1002> 

WRITE O.IOOS) 

DO 49 I-l.K 

DO 49 J ■ 1 » K 
49 MR I TE C 3 * 1 OO 1 ) R< t . J ) .X ( 1 • J > . I * J 

WRITEI3.1003) 

WRITE 1 3* 1006) 

DO 201 1-1 *K 

201 WRITE<3.1001 > TR( I » tTI ( I ) . I 

202 IFINCOOE.EQ.l )GO TO 55 
102 CALL EXIT 

100 WRITEO.1007) 
GO TO 202 

101 WRITE(3«10O4) 
DO 2000 I-l.K 
DO 2001 J"ltK 

2001 READ (1*1001) W.Z 
2000 READ < 1.1 001) W.Z 

RE AO (l.lOOO) NCODE 

GO TO 202 

ENO 



FORTRAN LISTING l*l0-F0-970 
SUBROUTINE MAT INV 
DIMENSIONAJ 25,26 ) ,B« 25,26) 
COMMON* »B,H 

01007 FORMAT ( 42H AN ELEMENT ON THE MAIN DIAGONAL IS ZERO ) 

KX*M*1 

00331*1. M 

0048L-1.M 

A(l,KX)=0.0 
00048 B(L,KXI*0.0 

A( I,KX)«l.U 

IV«1 

1FIAU, 1>.NE.0.0)G0T034 
IF(B( I,U.fcQ.O.O)GOTOlOO 

00034 Tl»A(I,l)»A(l,l>*B(I,l)»BU,l) 
TIR»A(I,1) 

T1X*B(I.I> 
0035J»1,KX 

WR»(A(I,J)»TlR*B(l. J)*T1X»/Tl 
Wl*IBU.J»»TlR-An.J)»TlXJ/Tl 
A( I. J) = MR 

00035 B(I,J)*WI 
IX*0 
IY-2 

IF( I.NE.MIG0T037 

00038 MX«I-l 
MY*1 
G0T039 

00037 MY=I*1 
HX-M 

00039 0040L*MY,MX 
IX*IX*l 
T1*A(L»1> 
T2-BIL. I) 
0040J-1.KX 

WR-A(I.J)*T1-B([.JI*T2 
WI«AII,J»»T2*BU.J1«T1 
A(L,JI»A(L, Jl-WR 

00040 B(L, J )*B< L.J )**MI 
IY«3 

IF « I.EQ.l )G0T044 
IFI(M-1).NE.IX)G0T038 
00044 D046LM.M 
0046 J- l.M 
NU»J*l 

AIL. J)«A(L,NUI 
00046 B(L,J)-BU,NU» 
00033 CONTINUE 

RETURN 
00100 HRITEI3.1007) 

STOP 

ENO 



APPENDIX VI 



FORTRAN LISTING 1410-FO-970 
DINENSIONRI 25,26), X (25,26), ERI25),E1( 25 ),TR(25),TI 125) 
COHNONR.X.K 

01000 FORMAT! IS I 

01001 F0RHAT(2F10.4,2I5) 

01002 FORMAT I /46H THE INVERTED MATRIX ELEMENTS ARE AS FOLLOWS! 

01003 F0RMAT!///46H THE SIMULTANEOUS EUUATION SOLUTIONS FOLLOW » 

01004 FORMATI 67H THIS MATRIX IS GREATER THAN 25X25. THE NEXT CASE WILL 

01005 FORMAT OOH REAL I MAG ROM COL) 

01006 FORMAT! 2 7H REAL I MAG ELEMENT I 

01007 FORMAT (42H AN ELEMENT ON THE MAIN 01 AGONAL IS ZERO) 

01008 FORMATI IX, S2H ) 
03000 FORMATI IH1 I 

0005* REAOIltlOOt) 

C READ IN THE SIZE OF THE MATRIX 

READ(1,1000)K 

IF(K.GT.25)G0T010l 

00501*1. K 

D030J-1.K 
00030 READI1,1001)R(I,J),XII,J) 
00050 REA0(ltl001)ER(t),EIII) 

C READ IN THE LAST COOE CARD. I FOR MORE CASES TO FOLLOW, AND 9999 

C FOR THE LAST CASE 

REA0!1,1000)NC0DE 

CALLMATINV 

D053IM.K 

TRI I )*0.0 

Till 1*0.0 

D053J«l,K 

Tl«R!I,J)*ERU)-X!I,J)*EHJ) 
T2*X! I,J)*ERI J)*RI I.J)*EI( J) 
TR(I)*TRII)+T1 
00053 Tim*TIII)*T2 
WRITEI3, 10001 
WRITEI3.1008) 
WRITE<3,1002) 
WRITfc(3,1005) 
00*91-1, K 
0049 J* I, K 

00049 WRITE 13,1001 )R( It J) ,Xl I , J ) , I . J 
WRITE! 3,1003) 
WR1TEI3.1006) 
002011-1, K 

00201 WRITEI3.100l)TR< I).TIIIJ,I 

00202 IFINCOOE.E0.DGOT055 
00102 CALLEXIT 

00100 WRITEI3, 10071 
GOT0202 

00101 WRITEI3, 10041 
G0T0102 

END 



© 



OUTPUT - CASK 1 



THIS IS A 4X4 MATRIX 

THE INVERTED MATRIX ELEMENTS AWE AS FOLLOWS 



SOLVED. 



EAL 


I MAG 


ROW 


COL 


- . 1 862 


1 .0951 


1 


1 


-.0085 


.0151 


1 


2 


-.0343 


.1282 


1 


3 


• 2154 


-1.2256 


1 


4 


.2296 


-1 .2040 


? 


1 


.0297 


-.0935 


2 


2 


-.O062 


.0331 


2 


3 


-.1673 


1.1311 


2 


4 


.2177 


-1.1513 


3 


1 


-.0080 


-.0049 


3 


2 


-.0056 


.0339 


3 


3 


-.1386 


1 .0672 


3 


4 


.0327 


-.0725 


4 


1 


.0197 


-.0779 


4 


2 


.0054 


-.0046 


4 


3 


.01 75 


.0319 


4 


4 



THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW 
REAL I MAG ELEMENT 

-.0137 .0129 I 

.0858 -.1333 2 

.0654 -.0551 3 

.0754 -.1231 4 



OUTPUT - CASE 2 



THIS I S A 3X3 MATRIX 

THE INVtRTED MATRIX rL E MP. N T S ARt AS FOLLOWS 



E*L 


I MAG 


ROW 


COL 


1 .841 7 


-5.3149 


1 


1 


.0385 


-.0488 


1 


'?. 


.0096 


-.0119 


1 


3 


.0385 


-.0488 


2 


1 


t .8432 


-5.3153 


2 


2 


.0385 


-.0488 


2 


3 


.0096 


-.0119 


3 


1 


.0385 


-.0488 


3 


2 


1 .841 7 


-5.3149 


3 


3 



THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW 
REAL I MAG ELEMENT 

•0748 -.3108 1 

-.2049 .5882 2 

-.2063 .5904 3 



c 

TECH 




1 026 2034 8 



